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We investigate the efficiency of systems of molecular motors operating at maximum power. We 
consider two models of kinesin motors on a microtubule: for both the simplified and the detailed 
model, we find that the many-body exclusion effect enhances the efficiency at maximum power of 
the many-motor system, with respect to the single motor case. Remarkably, we find that this effect 
occurs in a limited region of the system parameters, compatible with the biologically relevant range. 
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Molecular motors are cyclic machines that convert 
chemical energy into useful work, and perform a vari- 
ety of cellular functions such as transport, translocation 
and polymerization [l[ . 

Such machines work in environments at constant tem- 
perature, differently from, e.g., heat engines, whose effi- 
ciency is bounded by the Carnot's law. Thus, the effi- 
ciency of molecular motors is constrained by the thermo- 
dynamic limit 1, which can, however, only be achieved in 
the limit of vanishing power output. Driven by this obser- 
vation, several researchers have investigated the efficiency 
of microscopic engines working at maximum power, with 
the aim of characterizing the optimal balance between 
these two thermodynamic quantities While these 

works were focused on the efficiency of single motors, 
little is know on the efficiency of systems of interacting 
molecular motors. Indeed, molecular devices such as ki- 
nesin and myosin move on crowded filaments, where they 
can encounter other motors, and the resulting interaction 
can become non-negligible. In fact, molecular motor traf- 
fic is a relevant and widely studied phenomenon, which 
is typically modelled by using exclusion processes on lat- 
tices [i^oi. 

In the present letter we investigate the effect of motor- 
motor interaction on the system efficiency in the maxi- 
mum power regime, by using two different motor mod- 
els for kinesin under external mechanical force. The first 
model is characterized by a single internal conformational 
state: it is a reduced version of the lattice models that 
can be found in literature [l6j|, and it is mainly used to 
discuss in simple terms the issue of EMP in many-body 
systems. However, as emphasized in several works |9|- 



12|, it is crucial to include internal states if one wants 
to characterize properly the mechano-chemical cycles of 
molecular devices. Thus, we exploit the 6-state model for 
single kinesin introduced in [17j . and extend it in order 
to investigate the exclusion interaction effect on several 
motors moving on a lattice at the same time. 

For both the models we find the striking effect that the 
EMP is enhanced, as compared to the single motor case, 
by the mutual exclusion interaction. The increase in the 
EMP is mainly driven by the dynamical phase transi- 



tion that the system undergoes as the mechanical force 
is increased. Most remarkably, the increase of the EMP 
occurs in a narrow range of parameters, corresponding to 
the biologically significant one. 

Model I: Model with no internal states. The simplest 
description of a system of interacting molecular machines 
moving on a filamentous track is obtained by neglecting 
the internal conformational changes of the motor, and 
modelling the stepping of an individual motor by a single 
Poissonian step. The stepping is subject to an exclusion 
rule similar to the widely studied asymmetric simple ex- 
clusion process (ASEP) [15[. The motor is thus described 
as a particle moving on a lattice with a lattice constant 
a corresponding to the step size of the motor. The for- 
ward and backward jumping rates are denoted by p and 
q, respectively. In the absence of interactions the motor 
will thus exhibit a steady-state velocity vq = a(p — q). 
However, when the exclusion rule is introduced, the step 
is rejected, if the motor attempts to step into a site al- 
ready occupied by another motor, and the steady-state 
velocity v of a single motor is smaller than vq. The par- 
ticles bind to the filament at the left end with rate a and 
leave the track at the right end with rate /3. The mo- 
tor performs work against the external load force / < 
while hydrolyzing one ATP molecule into one ADP and 
one Pi molecule during each step. For constant concen- 
trations of ADP and Pi, we can write the transition rates 
according to chemical kinetics as 



P = w e 



(A/j+/a9)/T 



q = w e 



-fa(l-6)/T 



(1) 



where ujq is a microscopic rate, 9 is the load sharing 
factor, expressing the coupling between the force and 
the system kinetics, and T is the temperature, which is 
taken to be T — 4.1 pN • nm. Here and in the follow- 
ing we take ks = E The quantity Afj, is the chem- 
ical free energy of ATP hydrolysis, which can be ex- 
pressed in terms of the reactant concentrations: A/j, = 
Tln{A' oq [ATP]/([ADP][Pi])}, where K oq is the equilib- 
rium constant of the reaction. Thus, in this simplified 
model motor, the mechanical and the chemical cycles are 
tightly coupled, and an ATP molecule is hydrolysed only 
when a mechanical step takes place. The velocity of a 
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single motor is equal to v = J j p, where J is the steady- 
state probability current of the particles on the lattice, 
and p is the (average) bulk density of the motors. The 
quantities J and p can be obtained exactly, and the phase 
diagram of the probability current in the thermodynamic 
limit is that of a standard ASEP see also App. Al 
in [l~8[. When projected onto the effective parameters 
a = a/(p — q) and f3 = f3 / (p — q), it consists of three 
regions termed low-density (LD), high-density (HD) and 
maximal current (MC) phase, respectively. The current 
is given by J = ap{\ — ap), where p depends on the 
phase, and its value can also be obtained within several 
mean- field approaches [HI, [l^. In the LD phase, charac- 
terized by a < min(/3, 1/2), the density is p = a/ a, while 
for (3 < min(o!, 1/2) the system is in the HD phase and 
p = (1 — /3)/a. Finally, in the MC phase, characterized 
by a > 1/2 and /3 > 1/2, the bulk behaviour becomes in- 
dependent of boundary conditions and p = l/2a. Hence, 
the velocity of the interacting motors becomes 




aa for a < \, (3 > a (LD) 
for ft < i, P < a (HD) 
for a > \, > \ (MC). 



(2) 



It is worth noting that while the total particle current J 
is maximized in the MC phase, the single particle veloc- 
ity v attains its maximum in the LD phase as expected, 
since the hindrance due to the other particles is minimal 
in this phase. Moreover, the velocities fulfill the relation 
vo > «ld > "mc > "hd- In the HD phase the veloc- 
ity is constant due to queuing and only depends on the 
detachment rate p. Finally, we note that the typical mi- 
crotubule lengths L ~ 10/im [h| are long compared to 
the motor step size a = 8 nm, and the thermodynamic 
limit is thus applicable. 

Efficiency at maximum power. We proceed by evalu- 
ating the EMP for model I. We fix the value of Ap (which 
corresponds to fixing the ATP concentration as discussed 
above), and look for the optimal load force /* for which 
the output power P ou t = — fv is maximum, i.e. /* is 
solution of dP out /df = — (v + fdv/df) = 0, where v is 
determined from eq. ©. It is worth noting that chang- 
ing / changes the values of the forward and backward 
rates p and q, as given by eq. ([T]), and thus the values 
of vq as well as a and p. This in turn may induce a 
phase transition, according to the rules in eq. ([2]) , as will 
be discussed below. Since the motor is operating in the 
tightly coupled regime, the efficiency is simply given by 
r] = -Pout/ Pin = - fa/ A/i = f/f s , where P in = Apv/a 
is the input power, and f s = — Ap/a is the stall force 
Therefore, the EMP for a given value of A/i (and thus of 
f s ) reads 77* = f* / f s . We then repeat this procedure for 
increasing values of Ap. We choose the values of the four 
parameters cjo, Q, a and P as reported in fig. [1] caption 
Such values are suitable for the kinesin motor, see App 
A2 UM for a detailed discussion. 
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Figure 1. Model I: a) EMP for a single-state model of kinesin 
interacting through self-exclusion (solid lines) for two differ- 
ent values of the load distribution factor 9. For comparison, 
the EMP for non-interacting motors with the same parameter 
values is shown with dashed lines. Parameter values: a — 8 
nm, ujo = 1.3 ■ 10 -7 s _1 , a = 5 s~ , /3 = 3 s _1 . The values 
for 9 axe taken from [TtJ and correspond to two independent 
experiments on the kinesin motor. b)-d) P ou t as a function 
of / for three different values of the chemical free energy Ap 
and for 9 — 0.3. The vertical dashed line in c)-d) represents 
the phase boundary between the HD and the MC phase. 



The resulting EMP is shown in fig. [T^,, where two 
different regimes for the behaviour of such a quantity 
are observed. When the free energy of hydrolysis is 
small, the motor operates at maximum power in the MC 
phase, see fig. [TJ3. Since the velocity in this regime is 
v = i'mc = Vq/2, the resulting EMP is identical to that 
obtained for the single motor in the absence of inter- 
action. In particular, we recover the well-known linear 
response value for the EMP, rf = 1/2, as A/i/T 
[2j-|4|]. By increasing Ap further, the system is in the HD 
phase at small /, but then the maximal P out is achieved 
in the MC phase, see fig. [TJ;. However, there is a critical 
value A/i c for the free energy, that is achieved when the 
velocity at maximum power in the MC phase is equal to 
the velocity in the HD phase, the latter quantity being 
independent of Ap and /, see eq. (TJJ). For the present 
choice of parameters one finds Ap c ~ 18.5 T, see App. 
18| . When Ap > Ap c , the maximum power is no 



A2 in 

longer obtained when the system is in the MC phase, but 
rather at the boundary between the two phases, and the 
load force maximizing P out thus reads /* = /hd-mc, see 
figure [TJi. Since the velocity of non-interacting motors, 
vo(f), decreases with the force for all values of / < 0, 
one finds that the optimal force /g is smaller than /* 
when Ap > Ap c . Thus, for Ap > Ap c , we observe an 
increase in the EMP, 77* = /* / f Sl as compared to the non- 
interacting system, due to a change in the characteristic 
force- velocity relation v(f). Furthermore, the transition 
force /hd-mc is such that /hd-mc ->■ fs as Ap/T goes to 
infinity. Hence, 77* = /*//« ~* 1 f° r Ap/T -> 00 entirely 
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due to the self-interaction of the motors on the lattice. 
We note, however, that the limit rf — > 1 for Ap/T — > oo 
is unphysical, since in real systems the velocity does not 
grow unbounded with increasing Ap because of the dis- 
sipation. One would rather expect the EMP to decrease 
for very large Ap, as we will see for the more realistic 
model II. However, large values for the EMP, 77* > 0.6, 
are achieved for Ap moderately larger than Ap Cl before 
the saturation begins to take place. 

We find that the enhancement of the efficiency at maxi- 
mum power, as compared to the EMP for non-interacting 
motors, is primarily determined by the load factor 9, as 
shown in fig. UK, and as discussed in App. A2 in [3]. It is 
interesting to note, that for the value of the chemical free 
energy under physiological conditions, Ap/T ps 20 — 25 
[lij , the model predicts an increase in the EMP up to a 
factor of four. 

Model II: Model with internal states. We now con- 
sider the 6-state network model for single kincsin motors 
introduced in ref. [I?}, and depicted in fig. [31 and extend 
the model to describe traffic of kincsin motors interact- 
ing through steric interactions as above. The exclusion 
rule affects the transitions between chemical states 1 and 
4 since these transitions represent mechanical steps be- 
tween neighbouring sites. As in model I, the motors can 
attach to the microtubule at the left boundary with rate 
a and detach at the right boundary with rate (3. The dy- 
namics of the system can be obtained analytically within 
a mean-field approximation by employing the maximal 
current principle (MCP) [Io|. The MCP states that the 
mechanical probability current through the open system 
in the thermodynamic limit, J m , is given by [201 ] 

[ max Jf BC (p) for Pl > p r 
I min J^bc(p) for Pl < Pr, 

\P£[Pl,Pr] 

where pi and p r are the densities of the left and right 
reservoirs, respectively. The mechanical current in the 
bulk, «/™ B c(p), is obtained by solving the mean- field mas- 
ter equations for the corresponding homogeneous system 
with periodic boundary conditions (PBC), and homoge- 
neous particle density p, see eq. (14) in App. B2 [18j . 
We find that J™ B g has a single maximum, and accord- 
ing to the MCP the topology of the phase diagram is 
therefore identical to that of a standard ASEP with the 
three phases as in model I, as can be seen in fig. 2 in 
App. B2 [3]. The reservoir densities pi <r as a function 
of the transition rates and the injection and removal 
rates a and /3 are obtained by equating the bulk current 
J™ BC with the particle currents at the boundaries, i.e. 
Jp BC (p;) = a(l — api), and Jp^c(Pr) = I3ap r . Here we 
have assumed that the motor can detach from the track 
from any internal state. 

The maximal current principle yields the bulk densi- 
ties p = pi in the LD phase, p = p r in the HD phase, and 
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Figure 2. Model II: The 6-state network model from [TtJ for 
a single kinesin motor is adapted to describe kinesin traffic by 
introducing self-exclusion. States marked by D or T represent 
a head domain containing bound ADP or ATP, respectively, 
while an empty circle denotes an empty head domain. The 
motor uses ATP hydrolysis to perform a forward mechanical 
step from site j to site j + 1 through the F + = 1 12341 ) dicycle, 
or a backward step to site j — 1 through B + = |45614), if the 
respective neighbouring site is unoccupied. Furthermore, the 
network contains the futile hydrolysis dicycle 1 1234561). 



p = p* in the MC region, where p* = max p Jp BC (p) de- 
notes the density that maximizes the current under PBC. 
The mechanical and chemical probability currents are 
given by J m = Jp B c(p) anc ^ ^° = ^pbcC/ 3 )' respectively, 
where p is the appropriate bulk density for the phase, 
and J PBC is the chemical probability current through 
the homogeneous system. The phase boundary between 
the LD and MC phases given by a c is obtained by solv- 
ing pi(a c ) = p* , while the critical value /3 C defining the 
boundary between the HD and MC phases is determined 
by the relation p r {/3 c ) = p* . Finally, the LD-HD phase 
boundary is defined by the equation J^(a) = J^ D (P) 
that yields a linear relation, (3 = a/3 c /a c . For fixed a 
and /3, and for given values of Ap and /, we can thus cal- 
culate the phase boundaries analytically and determine 
the phase that the motor operates in under these condi- 
tions, and hence J m and J c . One thus obtains the motor 
velocity v = J m / p and the hydrolysis rate r = J c /ap. 

Efficiency at maximum power. Since the motor 
is loosely coupled due to the presence of backward 
hydrolysis-driven stepping and futile hydrolysis, the ef- 
ficiency is now given by 77 = —fv/Apr, and generally 
f s < A/x/a. The efficiency at maximum power predicted 
for model II is shown in fig. [3^ for two different sets of 
parameters. The difference between these two sets con- 
sists only in different values for the chemical load factors 
X1.2, which express the coupling between the external 
mechanical force /, and the chemical transition rates, 
see ref. fl7| and App. Bl in [18]. For the parameters ob- 
tained in [17[ (squares in fig. by fitting the Visscher 
et al. 's experimental data [2l|] , we observe that the EMP 
is enhanced significantly, as compared to the EMP for 
non-interacting motors under the same conditions, when 
the input free energy is in the range Ap/T « 18 — 21 [22I ]. 
The mechanism for the EMP enhancement is well exem- 
plified by fig. [SJd where we compare P ou t for interacting 
and non-interacting motors. At low force the interacting 
system is in the HD phase, while at intermediate force 
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Figure 3. Model II: a) The EMP for interacting (solid lines) 
and non-interacting (dashed lines) kinesin motors as a func- 
tion of the hydrolysis free energy A/x for two different parame- 
ter sets. Squares: parameters obtained in [l7| for the Visscher 
et al. experiment [2l[. Crosses: as before but with chemical 
load parameters Xi,2 as in [23|] (see text). For all curves we 
take the biological concentrations [ADP] = 0.07 mM, [Pi] = 1 
mM [24| and vary the ATP concentration. Where relevant, 
we take a = 5 s _1 and /3 = 3 s~ . The biologically relevant 
regime A/x/T « 20 — 25 is marked with vertical lines. The dis- 
continuity in the EMP for interacting motors at A/x ~ 14. 5T 
arises due to a discontinuity in the maximizing force /* as a 
result of the MC-LD phase transition, b) P out as a function 
of / for the second parameter set and for A/x = 20T. As 
/ increases, the motor goes from operating in the HD phase 
to the LD phase, while the MC phase is entered close to stall 
force. The maximum power is achieved at the HD-LD bound- 
ary, c) The efficiency n(f) corresponding to the parameters 
as in b). The shift of /* to higher values results in a boost of 
the EMP, 77* , with respect to the non-interacting case, tjq. 



we observe a phase transition to the LD phase, and our 
model molecular machines operate at maximum power at 
the transition between the two phases. Thus the force /* 
maximizing T out is shifted to higher values with respect 
to the non-interacting system values /q because of the 
change in the force- velocity relation v(f). Since the effi- 
ciency, Tl{f)i m the LD phase is approximately equal to 
the efficiency for the non- interacting motors, rjo(f), the 
higher value for /* implies an increase in the efficiency 
at maximum power as long as A/x < 21T, as shown in 
fig. |3fc. When A/x f» 21T, the optimal force /* becomes 
so large that rf = rj(f*) = r?o(/o)- For even higher values 
of A//, the value of /* approaches the stall force, while 
the input energy is dissipated through the futile cycles, 
and hence if —> as A/x — > 00. 

The range of values for the free energy A/x for which the 
EMP exhibits an enhancement is seen to partly overlap 
with the biologically applicable regime A/x/T ~ 20 — 25 
[H]. We find that the width of the "peak" that deter- 
mines the region for the increase in the EMP, is very 
sensitive to the values of the chemical load parameters 
Xi,2) and not to variations of the other model parame- 
ters as discussed in App. B3 jl^|. As a second set of 



parameters we take the values Xi — 0.3 and X2 = 0.4, 
as obtained in 1231 b y fitting the data from the Carter et 
al.'s experiment |25| to a 7-state model that is an exten- 
sion of the 6-state model considered in the present work. 
We plot in fig. [3^, the resulting EMP for these values, 
while keeping all the other parameters fixed (crosses) . In 
this case, the EMP exhibits approximately a two-fold in- 
crease in the presence of interactions for a wide range 
of free energies, A/x/T ~ 18 — 40, with a wider overlap 
with the biologically relevant regime A/x/T ~ 20 — 25. 
We note that the values for the parameters Xi exhibit a 
strong dependence on the underlying model used to fit 
the experimental data (T7II23I]. Furthermore, the param- 
eters Xi dep end strongly on the experimental conditions: 
in ref. [17j, by fitting the same model to two different 
experimental data sets, [2l| and [25|, xi takes the values 
0.05 and 0.15, respectively. We note that other factors 
may also contribute to shifting the peak of EMP to higher 
values of A/i as discussed in App. B3 [l8| . 

Conclusions. We have studied the EMP for two dif- 
ferent models of kinesin traffic, and showed that the 
motor-motor interaction plays a fundamental thermody- 
namic role by enhancing the system efficiency at maxi- 
mum power, as compared to the single motor case. Our 
findings are amenable of experimental verification, given 
that a few experimental groups have started measuring 
the dynamical properties of many- motor systems [26, 27 1, 
which can be easily compared to the single-motor ones. 
We note that when multiple motors are involved in trans- 
porting a single cargo, the force- velocity characteristic is 
modified, as compared to the single-motor-single-cargo 
system @, HI, 28 1. However, in a recent in vivo study 
[29|, Leidel et al. conclude that the motors share the 
load force equally. Based on this, one can argue that 
groups of motors bound to several fluid cargos may ex- 
hibit a dynamical phase transition as the one discussed 
above, leading to a similar enhancement of the EMP. Fi- 
nally, we expect our findings to hold beyond the case 
of molecular motors on a lattice: in a generic system of 
mechanochemical motors, interacting with a general po- 
tential, one might presumably observe an enhancement 
in the EMP, provided that the system exhibits a dynam- 
ical phase transition between two regimes with different 
characteristic response to the external mechanical driv- 
ing. This is worthy of future investigation. 
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APPENDIX A - MODEL I 
Al. The phase diagram of the probability current. 

In the thermodynamic limit, N — > oo, the exact phase 
diagram of the current J consists of three regions when 
projected onto the effective parameters a = a/ (p—q) and 
j3 = /3/(p — q), see fig. 1. The phases are characterized 
as follows: 

• Low- density phase (LD): In this regime character- 
ized by a < min(/3,l/2), the injection process is 
rate-limiting. The bulk density is p = a/a and the 
average current is J = (p — q)a(l — a). 

• High-density phase (HD): When /? < min(<5, 1/2), 
the density p = (1 — /3)/a and the current J = 
(p — <?)/?( 1 — (3) are functions of the rate-limiting 
extraction step. 

• Maximal current phase (MC): For a > 1/2 and 
[3 > 1/2, the bulk behaviour is independent of the 
boundary conditions, and one finds p = l/2a and 
the maximum possible current J — {p — g)/4. 

The low- and high-density phases are separated by the 
line a = j3, across which the bulk density is discontin- 
uous. The density profile on this line is a mixed-state 
of shock profiles interpolating between the lower density 
p = a/ a and the higher density p — (1 — /3)/a. Both 
the bulk density and the current are, on the other hand, 
continuous across the transitions from the low- and high- 
density phases to the maximal current phase. Finally, 
it should be stressed that the uniform densities stated 
above are only asymptotically accurate in the bulk, but 
must necessarily vary across a narrow region near the 
boundaries of the lattice in order to satisfy the boundary 
conditions imposed by the extraction and injection in the 
steady-state. Combining the above relations for J and p, 
and keeping in mind that vq = a(j> — q), one obtains the 
expression for the velocity v = J/p of the motors, eq. (2) 
in the main text. 

The problem of the ASEP can be addressed using sev- 
eral mean-field approaches such as recursion relations, 
hydrodynamic equations or the maximal current hy- 
pothesis [1^4]. By introducing the dynamical variable 
rii = 0, 1, expressing the occupation of the i-th site, the 
average density is given by pi = (n,i)/a. Common to 
all the mean field approaches is the assumption that the 
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Figure 1. The phase diagram of the probability current for 
model I is that of a standard ASEP in the thermodynamic 
limit. 



correlations between neighbouring sites can be neglected, 

(nin i ± 1 ) ps {m)(n l±1 ) = p l p i ± 1 . (1) 

The resulting mean-field master equations governing the 
evolution of the density pi arc 



d Pl (t) 

dt 
dpi ft) 

dt 
d/Qjyft) 

dt 



= Ji+x/2 ~ Ji-i/2 fori = 2,...,JV-l (2) 
= a(l - pi) - Ji+i/2 (3) 

= Jn-1/2 ~ PPN; (4) 



where the expression for the mean-field probability cur- 
rent Ji+i/2 is 



■h+i/2 = m(i - Pi+i) - - Pi)- 



(5) 



In the steady-state limit, pi = 0, the probability current 
J MF is independent of the site: J MF = = Ji-i/2- 

We note that the phase diagram of the probability cur- 
rent predicted by mean-field treatments of the problem 
happens to coincide with the exact phase diagram dis- 
cussed above, where the exact treatment requires to find 
a steady state solution for the occupation probability 
P(n 1 ,n 2 , ■ ■ .,n N ) [1, 2, 4]. 



A2. Dependence of the EMP on model parameters. 

Model I contains four parameters, namely the micro- 
scopic transition rate wo, the load factor 9, and the at- 
tachment and detachment rates, a and /3, respectively. 



2 



The values of a and are taken to be 5 s -1 [5] and 3 
s -1 [6], respectively. The value of 8 is sensitive to the 
choice of model and experimental setup. In the present 
case we use the values for 8 obtained in [6] by fitting- 
model II to experimental data from two independent sin- 
gle molecule experiments on kinesin. The value of uiq 
is fixed in such a way that the single motor velocity at 
zero force is the same as the single motor velocity for 
the more detailed model II. More precisely, we consider 
the 6-state model discussed in [6] and which model II is 
based on, and plot the single particle velocity vq at / = 
as a function of the chemical free energy A/i for the Viss- 
cher et al. parameters [7]. The curve is then fitted by 
v o(f = 0) = au (e A,I ' T — 1) in order to determine ojq 
for model I. The best fit is obtained for luq = 1.3e — 7 
s -1 . We note, however, that the simple single state 
model I does not capture the saturation of the velocity 
for A/x > 30fcsT that occurs due to dissipation of energy 
through futile cycles. 

Whenever the probability of attachment is larger than 
the probability of detachment from the track, i.e. a > /3, 
as is the case for our system, the motor will not operate in 
the low-density (LD) phase for any values of A/i and the 
load force /. The efficiency at maximum power (EMP) 
obtained within this model is therefore independent of 
a, since the latter does not enter the expressions for the 
probability current and bulk density in the high-density 
(HD) and maximal current (MC) phases. Furthermore, 
the critical value for the free energy, A^t c , that separates 
the two regimes for the behaviour of the EMP discussed 
in the main text, is robust against variations in ujq, f3 
and 8. The quantitative dependence of A/j, c on the pa- 
rameters is discussed below. Variations in /3 of one order 
of magnitude induce a change in A/i c by approximately 
2T, while a corresponding change in 8 only shifts Afx c by 
~ 0.1T. Changing ujq by one order of magnitude gives 
rise to changes in A/i c of about ~ 2T. The value of 8 has 
a big impact on ?7*(A/x c ) as discussed in [8]. Moreover, 
the load factor determines how fast rj* — > 1 for A/i — > oo, 
with small 8 giving fast convergence. Under more gen- 
eral conditions 8 is not constant but depends on the load 
force / and decreases with / [8] , and we therefore expect 
that the increase in the EMP due to self-exclusion is even 
more pronounced than indicated by the results presented 
in figs, la and 3a in the main text. 

The value of Afj, c separating the HD and MC dynam- 
ical phases is given implicitly in the following way. The 
location of the maximum in the MC phase, f*(Afi), is 
obtained by solving 



dP, 



MC 



df 



df [ 
d 



■fv M c(f,Atx)) 
(- §/(p(A/i, /)-«(/))) =0, (6) 



while Afi c is then determined by equating the velocity 
at maximum power in the MC phase with the constant 



velocity in the HD phase, 

VMc(f*(AfJ>c), A/i c ) = v H d = a/3- (7) 

The dependence of A/i c on the parameters of the model 
is thus logarithmic due to the relationship between the 
thermodynamic and mechanical quantities as expressed 
in the parametrization of the jumping rates p and q, eq. 
(1) in the main text. 

APPENDIX B - MODEL II 

Bl. Definition of the chemical load factors xi,2- 

The transition rates u>ij arc equal to the number of 
transitions from state i to state j per unit time. These 
rates can be parametrized as 



= Wy )0 $ij(/) with $ y (0) = l, 



(8) 



where uiijfi are the transition rates in the absence of ex- 
ternal load force, i.e. for / = 0. The force factors $y 
describe the effect of mechanical strain on the enzymatic 
reactions taking place in the motor domains [6]. For the 
mechanical transitions |41) and 1 14) the force factors take 
the form 



$4i(/) = e 9/a/T and $ 14 (/) 



-{l-9)fa/T 



(9) 



where 8 is the (mechanical) load factor already intro- 
duced in model I. For the purely chemical transitions, 
Liepelt and Lipowsky use the parametrization [6] 



$<,■ =2/(l + exp[-x ij '/a/T]), 



(10) 



where Xij are the chemical load factors. They find that 
the parameters Xij take on only two different values, and 
in the main text we therefore use the notation xi = X12 = 
X23 = X45 = X56 and % 2 = X34 = Xie- We note that the 
parametrization in eq. (10) is not standard and takes on 
a different form in e.g. [9]. 



B2. The phase diagram of the probability current. 

Solving for the mean-field dynamics of the open system 
coupled to reservoirs using the maximal current hypoth- 
esis (MCH) requires knowledge of the solution for the 
closed system with periodic boundary conditions (PBC), 
as can be seen from eq. (3) in the main text. The density 
of particles in the chemical state i is denoted p l , where 
i = 1, . . . , 6, see fig. 2 in the main text. The mean-field 
master equations governing the evolution of the densities 
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p are 

p 1 =(w 4 ip 4 - W14P X )(1 - p) + (u-2ip 2 - W12P 1 ) 





4- ft,),.-* rfi 


— ^lep 1 ) 








p = 


/ 2 
- (w 2 ip 


- U12P ) + 


(W32P 


-W 23 P ) 


1— 1 


p 3 = 


- (w 3 2P 3 


- W23/3 2 ) + 


/ 4 


~ W 34 p 3 ) 


O 


p 4 = 


C 4 


- Wi4/0 1 )(1 


-P)" 


(w 43 p 4 - w 34 p 3 ) (11) 






+ (w 6 ip 6 










p 5 = 


- (a; 54 p 5 


- o; 45/ 9 4 ) + 


(^65P 6 


- U56P 5 ) 




p 6 = 


- (w 6 5P 6 




(W 6 1P 6 


- W16P 1 ), 





where Uij is the transition rate between chemical state 
i and j. We note that p 4 is independent of the position 
on the lattice, since all lattice sites are equivalent un- 
der PBC. The total density p = ^\ p % thus serves as a 
control parameter. The mechanical current under PBC, 
JpBc(p), can then be found as 



idP) = (^41P 4 - Wi 4 p 1 )(l - p), 



(12) 



while the corresponding chemical current under PBC, 
Jp BC , is obtained from [10] 

Jpbc(p) = ( W 34P 3 - w 43 p 4 ) + (w 4 5P 4 - w 54 p 5 ). (13) 

A plot of JpB C (p) and Jp BC {p) for a specific set of pa- 
rameters can be seen in fig. 2. 

The MCH has proved to be a powerful tool in ana- 
lyzing exclusion processes since it is dictated by rather 
general dynamical properties of driven diffusive systems 
[4] . The microscopic details of the system only enter in 
so far that they determine the functional form of the cur- 
rent JpBc(/°) ano - * ne reservoir densities pi_ r . Because the 
probability current j of a standard ASEP under PBC ex- 
hibits a single maximum, since j (p) = p(l — p), the MCH 
predicts that the topology of the phase diagram for a 
given exclusion process is identical to that of a standard 
ASEP (fig. 1), whenever Jpbc(p) nas a sm gl e maximum 
as in our model II (see fig. 2). However, while we find 
the three phases encountered for model I, the behaviour 
of the critical points determining the phase boundaries 
becomes more involved. For given values of Ap and /, we 
can nevertheless analytically calculate the phase bound- 
aries that depend on all the transition rates of the model 
and hence determine the phase that the motor operates 
in under these conditions, when the attachment and de- 
tachment rates from the track are fixed. An example of 
the resulting phase diagram when A/i = 20T and / = 
is provided in the inset of fig. 2. 



B3. Dependence of the EMP on model parameters. 

The symmetry of the forward and backward stepping 
cycles of the model implies that the transition rates uiij 




Figure 2. Model II: The mechanical and chemical probability 
currents, Jpbc an d Jpbc> as a function of the occupation 
probability ap for interacting kinesin motors on a lattice with 
PBC for Ap = 20T, / = and the parameters obtained in [6] 
for the Visscher et al.'s experiment [7]. Note that the currents 
are not symmetric around ap = 0.5 as for a standard ASEP 
due to the lack of particle-hole symmetry in model II. Inset: 
The phase diagram of the current obtained from Jpbc using 
the MCH. 



fulfill, 



W45 
^56 
W61 



W12 
W23 
W34 



W54 = W21 
^65 = W32 



w 16 — (^14,o/^41,o) 2 W43) 



where the relation for ujiq follows from the local detailed 
balance conditions, cq. (3) in [6]. This leaves us with 
eight independent rate constants, W12, OJ21, ^23, ^32-, ^34, 
W43, W41 and W14. In order to preserve the local detailed 
balance conditions, these rate constants were varied pair- 
wise, such that u>ij and u>ji were changed simultaneously 
by the same amount. In the following, all variations of 
the parameters refer to changes of one order of magni- 
tude. We study semi-quantitatively the effect of such 
variations on 1) the position of the "peak" where the 
EMP enhancement occurs, 2) the width of the peak, and 
3) the (maximum) magnitude of the amplification of the 
EMP as compared to the single-motor case, see fig. 3a 
in the main text. We will now discuss the robustness 
of these three measures to variations in the four transi- 
tion rate pairs, the load factor 9, the rates a and /3, and 
the concentrations [ADP] and [Pi]. The influence of the 
chemical load factors xi,2 is addressed in the main text. 

An increase (decrease) in the rates 6^34,43 shifts the po- 
sition of the peak shown in fig. 3a in the main text to the 
left (right) by about 2T, while leaving the width and the 
amplification essentially unchanged. A decrease in 0^2,21 
or a>23,32 only decreases the amount of amplification to 
approximately 10%, while an increase in the same pa- 
rameters does not alter the EMP significantly. The same 
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is true for the rates 0^14,41, besides that the amplification 
in this case is reduced to about 40%. A gradual increase 
in the load factor starting from 9 = 0.3 lowers the ampli- 
fication until it eventually disappears for 9 close to 0.7. 
Gradually decreasing 9 from 9 = 0.3, however, does not 
have a significant effect on the EMP. 

Increasing the attachment rate a slightly increases the 
amplification, while the enhancement of the EMP disap- 
pears altogether when decreasing a (or raising the de- 
tachment rate /3). The latter is due to the fact that the 
system is never in the HD phase whenever a < /?. How- 
ever, although some disagreement is present in literature 
with regard to the actual values of a and f3, most refer- 
ences suggest that a > (3 for kinesin motors (see [11, 12] 
and references therein). Decreasing (3 implies a shift of 
the peak to lower values of Afj, by approximately 2T, 
while increasing /3, with the constraint a > /3 satisfied, 
shifts the peak to higher values of Afi. A higher value for 
/3 as suggested in [13] would thus contribute to moving 
the region of EMP enhancement further into the biolog- 
ical parameter range. Furthermore, increasing (decreas- 
ing) the concentrations [ADP] and [Pi] moves the position 
of the peak to the left (right) by roughly 2T. The cellular 
concentrations of ADP and Pi vary over several orders of 
magnitude and are in general poorly known. Smaller val- 
ues for [ADP] and [P;] would decrease the transition rates 
of the futile hydrolysis and reverse stepping cycles of the 
motor, thereby resulting in a shift of the EMP curves to 
the right in fig. 3a in the main text. 
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Figure 3. Model II: The occupation probability ap as a func- 
tion of the load force / for three different values of the free en- 
ergy A/i within the physiological regime: Afi = 20T (crosses), 
23T (squares), 25T (dots). The vertical lines indicate the 
value of the force, /hd-ld(A/i), at which the phase transition 
between the HD and the LD phases takes place. Parameter 
values are as in fig. 2 but with xi = 0.3 and X2 = 0.4, and 
correspond to the second parameter set in fig. 3a in the main 
text. 



B4. Motor density under physiological conditions. 

When the condition of maximum power is relaxed 
down to approximately 90% of maximum power, the sys- 
tem will operate in the low-density phase under physi- 
ological conditions, A/i/T ~ 20 — 25 [14]. In this case, 
the motor density on the track will lie within the range 
0.4 - 0.5, sec fig. 3. 
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